MIT-CTP 4177 
INT-PUB- 10-042 

Radiation of a circulating quark in strongly 
coupled JV = 4 super Yang-Mills theory 

Christiana Athanasiou*, Paul M. Chester*, Hong Liu*, Dominik Nicked 

and Krishna Rajagopal* 

*Center for Theoretical Physics, MIT, Cambridge, MA 02139, USA 
^Institute for Nuclear Theory, University of Washington, Seattle, WA 98195, USA 

Abstract. The energy density and angular distribution of power radiated by a quark undergoing 
circular motion in strongly coupled jY = 4 supersymmetric Yang-Mills (SYM) theory is computed 
using gauge/gravity duality. The results are qualitatively similar to that of synchrotron radiation 
produced by an electron in circular motion in classical electrodynamics: At large velocities the 
quark emits radiation in a narrow beam along its velocity vector with a characteristic opening angle 
a ~ 1/y and radial thickness scaling like ~ 1/y 3 - 
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INTRODUCTION AND MOTIVATION 

Through the discovery of gauge/gravity duality [1], the dynamics of quantum fields in 
certain strongly coupled gauge theories can be addressed by the dynamics of classical 
fields living in a higher dimensional spacetime. In particular strongly coupled real-time 
phenomena in ,jV = 4 super Yang-Mills (SYM) theory have become accessible to reli- 
able calculation for the limit of large colors N c and large but finite t'Hooft coupling X. 
Much of the analysis in this direction has so far been focused on finite temperature sys- 
tems, whose strongly coupled plasmas are expected to behave similarly to the strongly 
coupled quark-gluon plasma of QCD being produced and studied in heavy ion collision 
experiments at the Relativistic Heavy Ion Collider [2-5]. For reasons that will be out- 
lined below, we will however consider the specific case of a test quark in circular motion 
in the vacuum of SYM theory, i.e. synchrotron radiation. For this case we are able to 
find an analytic result that can be confronted with expectations from the weakly coupled 
regime, electrodynamics and a parton picture. 

Because the vacuum of SYM theory is not similar to that of our world, studying 
the radiation of a test quark in circular motion in this theory does not have direct phe- 
nomenological motivation. However, we shall find that the angular distribution of the 
radiation is very similar to that of classic synchrotron radiation: when the quark is mov- 
ing along its circle with an increasingly relativistic velocity, the radiation is produced in 
an increasingly tightly collimated beam — a beam which, if Fourier analyzed, contains 
radiation at increasingly short wavelengths and high frequencies. 
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First, this finding is interesting in its own terms. From the point of view of the non- 
abelian gauge theory it is not obvious why the angular distribution of the radiation at 
strong coupling is similar to what is known for the weakly coupled regime. Indeed 
a study related to zero temperature radiation analyzed decays of off-shell bosons in 
strongly coupled holographic conformal field theories [6] and it was shown that there 
is no correlation between the boson's spin and the angular distribution of power radi- 
ated through a sphere at infinity. In other words, if one prepares a state containing an 
off-shell boson with definite spin, in the rest frame of the boson the event averaged an- 
gular distribution of power at infinity will always be isotropic. Similar behavior regard- 
ing the isotropization of radiation at strong coupling was also reported in Refs. [7, 8]. 
A reasonable question is whether the isotropization is a characteristic of the particular 
initial states studied in Ref. [6] or whether it ia a natural process which happens dur- 
ing the propagation of strongly coupled radiation. It has been suggested in Ref. [7] that 
the mechanism responsible for the isotropization is that parton branching is not sup- 
pressed at strong coupling and, correspondingly, successive branchings can scramble 
any initially preferred direction in the radiation as it propagates out to spatial infinity. 
If the intuition about the emission and propagation of radiation derived from these cal- 
culations applies to radiation from an accelerated test charge, our result should show 
isotropization of the energy density when propagating towards infinity. However, this is 
not found. 

Second, our results in this formal setting nevertheless open the way to new means of 
modeling jet quenching in heavy ion collisions. If we were to add a nonzero temperature 
to our calculation, we could watch the tightly collimated beam of synchrotron radiation 
interact with the strongly coupled plasma that would then be present. The beam of 
radiation should be slowed down from the speed of light to the speed of sound and 
should ultimately thermalize, and it would be possible to study how the length- and 
time-scales for these processes depend on the wavelength and frequency of the beam. 
The authors of Ref. [9] have shown that the power that it takes to move the quark in 
a circle is given by the vacuum result [10] even for a test quark moving through the 
strongly coupled plasma present at finite temperature T as long as cofi y 3 3> K 2 T 2 . Here 

a>o is the angular frequency, v the quark's speed and 7=1/ y/l — v 2 . In the opposite 
regime, where co^y 3 <C K 2 T 2 , the power it takes to move the quark in a circle through 
the plasma is the same as it would be to move the quark in a straight line. This power has 
been computed using gauge/gravity duality in Refs. [11, 12], and corresponds to pushing 
the quark against a drag force. The analogue of our calculation — namely following 
where the dissipated or "radiated" power goes — was done, again using gauge/gravity 
duality, in Refs. [13-16]. It was found that the radiation pattern eventually takes the 
form of hydrodynamic excitations at distances large compared to l/T. For an extension 
of the present calculation to nonzero temperature, we therefore expect synchrotron 
radiation for ft^y 3 3> n 2 T 2 at length scales that are small compared to l/T, which 
should eventually thermalize locally, likely converting into an outgoing hydrodynamic 
wave moving at the speed of sound, before ultimately dissipating and thermalizing 
completely. Watching these processes occur may give insight into jet quenching in the 
strongly coupled quark-gluon plasma of QCD. 
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FIGURE 1. Cartoon of the gravitational description of synchrotron radiation at strong coupling: The 5d 
geometry of AdSs is shown as a foliation with leaves isomorphic to Minkowski space and understood as 
slices in the figure. The boundary as an example, is mapped out with color-coding indicating the induced 
energy density at a given time. The string (purple line) is attached to the boundary corresponding to the 
location of the quark in the dual theory and hangs down into the bulk. 

GRAVITATIONAL SETUP 

The determination of the energy radiated by a moving test quark (specifically, an in- 
finitely massive spin 1 /2-particle of a JV = 2 hypermultiplet) coupled to SYM theory 
in the large N c limit for large but finite X can be translated into a classical gravitational 
description in an AdSs background using the gauge/gravity dictionary [1, 17]. A cartoon 
of the the dual realization is shown in Fig. 1: The four dimensional boundary of AdSs 
has the geometry of ordinary Minkowski space. Ending at the boundary is a classical 
string whose endpoint follows a trajectory that corresponds to the trajectory of the quark 
in the dual field theory. The specific condition in our setup that the endpoint rotates about 
a specific (Minkowski space) axis results in the entire string rotating about the same axis 
with the string coiling around on itself over and over. Since we consider a single test 
quark in the large N c limit, we can neglect the backreaction of the string profile onto 
the geometry. However, the rotating motion of the string creates gravitational radiation 
which propagates up to the boundary of the geometry. Just like electromagnetic fields 
induce surface currents on conductors, the gravitational disturbance near the boundary 
induces a 4 d stress tensor on the boundary [18, 19], which corresponds to the expecta- 
tion value of the stress tensor in the dual quantum field theory. So, by doing a classical 
gravitational calculation in five dimensions (which happens to be somewhat analogous 
to a classical electromagnetic calculation) we can compute the pattern of radiation in the 
boundary quantum field theory, including all quantum effects and working in a strong 
coupling regime in which quantum effects can be expected to dominate. As indicated in 
Fig. 1, the induced Ad stress tensor inherits the coiled structure of the rotating string. 



In the following we will use the convention 



ds 2 = -J [-dt 2 + dr 2 + du 2 ] , ( 1 ) 

n 

for the metric of the AdSs spacetime, where L is the AdS$ curvature radius, u is the 
(inverse) AdSs radial coordinate with the boundary located at u = 0. We first outline the 
determination of the string profile, which then serves as an input for the calculation of 
the induced energy pattern. For details we refer to Ref. [20]. 



THE ROTATING STRING 

The dynamics of classical strings are governed by the Nambu-Goto action 



/A 

5 ng = ~ttzp> / dxdo^/-g, (2) 



2nL 2 

where a and x are the coordinates on the worldsheet of the string, and g = detg a f, 
where g ab is the induced worldsheet metric. The string profile is given by an embedding 
X M (x,o) of the worldsheet into AdSs. As we are interested in quarks which rotate at 
constant frequency (Oq about the z axis, it is convenient to choose x = t and o = u. 
Furthermore the possible set of solutions can be constraint to the parameterization 

X M (t,u) = (t,r s (t,u),u), (3) 

where in spherical coordinates {r, 0, <p} the three- vector r$ is given by 

r s (t,u) = (#(w),f,<K") + «oO • ( 4 ) 
The Nambu-Goto action then simplifies to 

a 



5ng = — ^ — / dtduJ2? , (5) 
2k J 



where 



J(1-(0 2 R 2 )(1+R' 2 )+R 2 ^ 2 
^ = 1 ~ 2 • (6) 

To determine the shape of the string in Fig. 1, we must extremize Eq. (5) in R(u) and 
<j)(u). 

For real- valued profile parameterizations R(u) and <j>(u) , it turns out that the solutions 
for the corresponding Euler-Lagrange equations are uniquely given by [20] 

R(u) = \Jv 2 y 2 u 2 +R-1, 0(«) = -uycoo + arctan (uycoo) (7) 



where v = RqCOq is the velocity of the quark and 7=1/ \J\—v 2 . 



As a by-product the shape of the rotating string also determines the total power 
radiated, as we now explain. The string has an energy density 7Tq and energy flux tZq, 
where 

with Sng given in (2) and where energy conservation on the string worldsheet requires 
d a ^M = 0. The energy flux down the rotating string then turns out to be 



K M = g ,l -\rft/f\ (8) 



P = -4 = ^a\ (9) 



where a = vco^y 2 is the proper acceleration of the test quark. The energy flux down the 
string is the energy that must be supplied by the external agent that moves the test quark 
in a circle; by energy conservation, therefore, it is the same as the total power radiated 
by the quark and agrees with the findings in Refs. [9, 10]. 



GRAVITATIONAL PERTURBATION AND INDUCED ENERGY 

DENSITY 

The starting point for the determination of the induced stress-energy tensor are the 
linearized Einstein equations 

2 P ^ 

— D h MN + 2D D( M h N ) P - D M D N h + j^hMN 

+ (^D 2 h - D P D® hpQ -^h^j G$ N = 2k 2 t MN , (10) 

where the full metric Gmn = G^ N + Hmn is linearized around the unperturbed metric 

G^ N defined in Eq. (1), h = h M M , Dm is the covariant derivative with respect to G$ N , K 2 
is the 5d gravitational constant and tMN is the 5d stress tensor of the string. 
According to the gauge/gravity dictionary, the on-shell gravitational action 

Sg= 2k£ / d5x ^(^ + J?) + S gh (ID 

is the generating functional for the boundary stress tensor [17, 19]. Here, G is the 
determinant of Gmn, & is its Ricci scalar, and Sqh is the Gibbons-Hawking boundary 
term [21], discussed and evaluated in the present context in Ref. [16]. The source of the 
dual stress-energy tensor is identified with 

u 2 

gfiv(x) = lim -J G M v (•*,"), (12) 
such that the boundary stress tensor is then given by [19] 



(13) 



with g denoting the determinant of g^ v . The gravitational perturbations are sourced by 
the 5d string stress tensor, for which in general 



f MN = _J* V ^ g ab daX M dbX N 5 3 {r _ rs) {u) 

Focusing on the energy density <S"{t,r) = T 00 (t,r), it is then possible - though non- 
trivial - to derive [20] 

g{t,r) = — d A r' / du9(t- 1') Uutw - t M5 V M W) — 
% J Jo \- U L 

+|r - rf (4t 00 - At 55 + 2f«) — ^ - (t^WVjW) 

where 



3u v J 1 J ' 2u 



(15) 



W = -(t-t') 2 + u 2 + \r-r'\ 2 (16) 

and V' M is the partial derivative with respect to X' M = {t\v\u}. Because we have ana- 
lytic expressions for the profile of the rotating string, it turns out that we can furthermore 
evaluate £ explicitly, albeit in a tedious calculation for which we again refer to Ref. [20]. 
We then arrive at the central result of this work: 

<^> r ) = 7T A 7.4^6 4ry 2 S(? ret -0+7ry 2 CT 2 S(? ret -0 3 +4y 2 ftj 2 (? ret -0 4 (17) 



247T 2 y 4 r 6 S 6 



+ (27 2 -4r 2 v 2 7 2 (O6 2 sin 2 0+3r 2 7 4 GJ o 2 S 2 )(? ret -O 2 -2r 2 S 2 
+ SvfcOor(t Kt -t)(t rst -t+rZ)smGcos((p-coot rst )] , 



where 

|r-r„(,,0)|-r.r„('.0) 
|r| 

and 

t-t Tet -\r-r s (t TeU 0) \ =0. (19) 

We have been able to obtain an analytic expression for our result in terms of the retarded 
time t Tet . For the case of a charge in circular motion, t ret must be evaluated numerically, 
similar to the corresponding case in classical electrodynamics [22]. 

Our results simplify in the far zone |r| — > <*>. With t, = 1 — v sin sin(<p — C0Qt iet ), here 
we find 

i r i2 r) co 2 yi 4-4v 2 sin 2 6-7^+3^y 2 

and from the continuity equation f ■ S(t, r) = $ (/, r) for the asymptotic energy flux. This 
allows us to obtain the time-averaged angular distribution of power 

dP _ v 2 r 3 ^v / l5r 2 -l-v 2 7 2 cos(20) 

\Zl ) 



da 327T 2 ( 7 2 cos 2 + sin 2 )5/2- 




v = 3/4 





FIGURE 2. Left: A cutaway plot of r 2 S/P for v = 1 /2. Right: A cutaway plot of r 2 S/P for v = 3/4. 
In both plots the quark is at jc = Rq, y = at the time shown and its trajectory lies in the plane z = 0. The 
cutaways coincide with the planes z = 0, <p = and <p = 7n/5. At both velocities the energy radiated by 
the quark is concentrated along a spiral structure which propagates radially outwards at the speed of light. 
The spiral is localized about = n/2 with a characteristic width 80 ~ 1/y. As v — s- 1 the radial thickness 
A of the spirals rapidly decreases like A ~ 1 /y 3 . 



We observe that radiated power is focused about 6 = n/2 with a characteristic width 
~ 1/y and - upon integrating over all solid angles - we verify the total power radiated 
agrees with Eq. (9). 



ILLUSTRATION 

Fig. 2 shows two cutaway plots of the energy density <§ at strong coupling (multiplied 
by r 2 /P where P is the total power radiated) produced by a quark moving on a circle 
of radius Rq at velocities v = 1/2 and v = 3/4. The figure is obtained by evaluating 
(18) numerically. The motion of the quark is confined to the plane z = and at the time 
shown the quark is at x = Rq, y = 0, and is rotating counter clockwise. The cutaways 
in the plots show the energy density on the planes z = 0, (p = and (p = Ik/ 5 where 
<p is the azimuthal angle. As is evident from the figure, energy is radiated outwards in 
a spiral pattern as the quark accelerates along its trajectory. This radiation falls off like 
l/r 2 and hence has a constant amplitude in the figure and propagates radially outwards 
at the speed of light. The figure shows the location of the energy density at one time; 
as a function of time, the entire pattern of energy density rotates with constant angular 
frequency (Oq = v/Rq. This rotation of the pattern is equivalent to propagation of the 



radiation outwards at the speed of light. As seen by an observer far away from the quark, 
the radiation appears as a short pulse just like a rotating lighthouse beam does to a ship 
at sea. 

From Fig. 2 we also see that the outward going pulse of radiation that does not broaden 
as it propagates. This can be confirmed analytically. As time progresses, the pulses move 
outward at the speed of light, and an observer at large |r| sees repeating flashes of 
radiation. The figure underlines that even though the gauge theory is nonabelian and 
strongly coupled, the narrow pulses of radiation propagate outward without any hint of 
broadening or isotropization. We find the same result at much larger values of y, where 
the pulses become even narrower - their widths are proportional to I/7 3 , as for classical 
electrodynamics [22] and also as for weak coupling in SYM theory [20]. 

Not visible in Fig. 2 but certainly worth mentioning is the fact that the energy density 
becomes negative before and after the pulse for sufficiently large velocities. We attribute 
this to quantum effects, showing that the classical gravitational setup is indeed dual to a 
strongly coupled quantum field theory. 
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